Cosmic Neutrino Flavor Ratios with Broken — v T Symmetry 
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Reactor neutrino experiments have now observed a nonzero value for #13 at 5a, and global fits to 
data imply a nonzero value above 10u. Nonzero values for #13 and/or 632- j break a v^-v r symmetry, 
which has qualitative as well as quantitative implications for the time-evolution of neutrino flavors. 
In particular, the large-distance flavor evolution matrix, non-invertible with v^ L -v T symmetry, is now 
invertible. This means that measurements of neutrino flavor ratios at Earth can now be inverted to 
directly reveal the flavor ratios injected at cosmically distant sources. With the updated values of 
the three neutrino mixing angles, we obtain the inverted large-distance evolution matrix and use it 
to derive several phenomenological relations between the injection flavor ratios and the observable 
ratios at Earth. Taking the three popular injection models as examples, we also exhibit the shift of 
Earthly observed flavor ratios from the corresponding values in the case with v^-v T symmetry. 

PACS numbers: 14.60.Pq, 14.60.Lm, 95.85.Ry 



I. INTRODUCTION 



It is well-known that a statistical average over a neu- 
trino ensemble from cosmic distances eliminates the 
quantum-mechanical phase (frjk = L(m 2 — m 2 .)/2E 
between states, leaving a relatively simple result for 
neutrino- flavor evolution. The evolution v a —¥ vp, with 
a and /3 any elements of the three-flavor set {e, /i, t}, 
is described by the propagation matrix P whose positive 
definite elements are 



P a p = J2\U aj \ 2 \U Pj 



(1) 



The physics behind this formula is that the correct basis 
for particle propagation is the mass basis, because the 
particle propagator in field theory is an analytic function 
with poles at mass values; we sum over the unobserved 
mass states labeled by j, and we weight each such mass 
state by its classical probability |f/ Q:) | 2 to overlap with 
the flavor a produced at the source, times its classical 
probability \Upj\ 2 to overlap with the flavor /3 detected 
at Earth. The sum on j = 1, 2, 3 is over the three active 
neutrino states. (Astrophysical distances are so much 
larger than relevant oscillation lengths that subtleties in 
the definition of cosmic distance [l| are not important 
here.) 

Phase-averaging restores CP-invariance, and so the 
matrix P describes both neutrino and anti-neutrino prop- 
agation. Furthermore, CP-invariance, according to the 
CPT-theorem, implies also T-invariance, and so the ma- 
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trix P is symmetric, P a p = Pp a . Explicitly, one has 



/ Ei \Uejf \U ej \ 2 Wei? \u M \ 2 E, \u e3 \ 2 \u TJ \ 2 \ 



p = 



V 



w = PW . 



(3) 



If P is an invertible matrix (i.e., has a nonvanishing 
determinant), then the inverse equation 



w = p- 



w 



(4) 



allows one to input the neutrino flavor ratios observed at 
Earth to obtain the flavor ratios dynamically injected 
at the cosmic source. At present, the only observed 
sources of extra-terrestrial neutrinos are from the Sun 
and SN1987A. The hope is that neutrino telescopes, re- 
cently deployed 0, [1], or soon to be deployed [J], will 
begin to observe neutrinos from more distant sources. 
When an ensemble of neutrino events are collected, track- 
topologies will allow one to glean the ratios of flavors 
arriving at Earth Very recently, the IceCube exper- 
iment has announced what is likely the first observation 
of high-energy extra-galactic neutrinos, two showering 
events characteristic of v e 's (or ^ e 's, since the experi- 
ments cannot distinguish v from v), with energies ~PeV. 
It appears that the era of neutrino astrophysics is sud- 
denly upon us. 

On the terrestrial neutrino experimental front, neu- 
trino mixing data were consistent with v fl -v T symmetry, 
defined operationally as |?7 W -| = \U. 
The popular example of such a v^- 
the TriBimaximal (TBM) model [f 



T] 1 , until very recently. 
u T symmetric model is 
] , which has the follow- 



(2) 

Thus, the flavor ratio unit-vector injected at the source, 
W = (W e ,W M ,W T ) is measured at Earth to be w = 
(w e ,Wfj,,w T ), where 
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ing classical probability and flavor propagation matrices: 



\U Q 



and 
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(5) 



(6) 



With exact v jl -v r symmetry, the second and third rows 
of the |J7| 2 and P matrices are identical, by definition. 
Accordingly, the determinant of matrix P vanishes, and 
the matrix is not invertible. Thus, prior theoretical 
work, which attempted to relate flavor ratios at Earth to 
the ratios injected at the cosmic sources, guessed at the 
source ratios and and evolved the guesses forward with 
Eq. ([3]) to obtain the observable ratios at Earthly detec- 
tors. Popular guesses are the pion decay-chain flavor- 
ratios |(1:2:0), the ,3-beam ratios (1:0:0) Q, and the 
incomplete pion decay-chain or quenched /x-decay ratios 
(0:1:0) @. We will examine flavor evolution in these 
three injection models in some detail below. A thought- 
ful overview of neutrino injection models is given in [9j . 

More satisfying would be to approach the study with 
observed flavor ratios, and evolve them backwards via 
Eq. (QJ to obtain as directly as possible the astrophysical 
quantities of interest, namely the flavor ratios injected at 
the sources. Recent neutrino data (Toj | from nuclear reac- 
tors reveals that the v^-Vr symmetry is broken. Hence, 
the determinant of P is no longer vanishing, and the in- 
verse flavor propagation matrix is calculable. 

The first purpose of this paper is to provide this inverse 
propagation matrix. As the second purpose of this inves- 
tigation, we will draw various phenomenological infer- 
ences for three- flavor neutrino astrophysics. For example, 
we will plot the movement of Earthly flavor ratios away 
from their TriBiMaximal values, for the three most popu- 
lar cosmic-source flavor models; present a relation among 
flavor ratios at Earth that determines whether tau neu- 
trino's are injected at the source; derive a general formula 
for the injection flavor ratio at the source in terms of the 
observable ratio of readily measurable track and shower 
events at Earth; and derive bounds on the possible fla- 
vor ratios to be observed on Earth, as implied by the 
measured mixing angles. If observations at Earth were 
to violate these latter bounds, then some physics other 
than flavor mixing via phase-averaged vacuum oscilla- 
tions would be at play. Examples of new physics could 
be neutrino decay [llj , or oscillations into new states such 
as sterile neutrinos or pseudo-Dirac states fl2l |. 



II. TWO-COMPONENT FLAVOR 
EXAMPLE 



It is instructive to see how flavor bounds may be de- 
rived in the simple case of TBM mixing. The v^-v T sym- 



metry tells us that and v T will arrive in equal num- 
bers, regardless of the flavor distribution at the source. 
Thus, there are but two relevant flavor ratios at Earth, 
w e and wy, with uy equally split between and v T , and 
w e + uy = 1. We may obtain the propagation matrix in 
this (e, /£) basis by adding the identical ty and v T rows in 
Eq. ©, and omitting the now redundant third column. 
One gets 



peff 
r TBM 



5 2 
4 7 



The propagation equation 



peff [ W e 



(7) 



(8) 



is linear, so the flavor extremes at Earth are found by 
inputting into Eq. J5]) the pure flavor vectors W — (1,0) 
and W — (0, 1) at the source. The results are trivially 
2/9 < w e < 5/9, correlated with 7/9 > uy > 4/9. Such 
is the allowed region of flavor space at Earth (w) in the 
TBM model. 

Proceeding further, the determinant of P^bm ^ s 
nonzero, and so this matrix is invertible. The inverse 
matrix is 



(P 



TBM) 



7 

-4 



We 



TBM/ 



uy 
(9) 

From Eq. © we derive an interesting expression relating 
flavor ratios at the source to the same ratio observed at 
Earth: 



We 

Wy 



7 



5-4 a* 



(10) 



Neutrino telescopes are particularly adept at distinguish- 
ing muon tracks due to interactions, from showering 
events due to v e and v T interactions. We may bin to- 
gether the latter events as w s h = w e + w T . Then, in- 
putting the TBM relations w e = w s h — uy and = 2uy 
into Eq. (fTU)) . we get an alternative expression of the same 
relation, 



We 

Wy 



7-11 



14 



(11) 



These TBM relations, Egs-lfTU]) and ([TT]), hold for any 
injection model. If either LHS were output, its value 
would discriminate among injection models. 

Theory strongly suggests that v r production at the 
source is very small |13l |14| (i.e., W T is zero or nearly 
THE TBM S0; or equivalently, that We> « W^). In hand with this 
assumption, the inverse matrix (Ptbm) -1 of Eq. (|9]) pro- 
vides a complete reconstruction of the flavor ratios at the 
source in terms of flavor ratios observed at Earth, for the 
TBM-mixing model. 



3 



III. BROKEN v^-Vr SYMMETRY AND 
THREE-COMPONENT FLAVOR 

With the observation of #13 far from zero, we learn that 
v^-v r symmetry is likely broken, and an analysis of the 
full three-component vectors w and W of Eqs. © and 
(|4"|) is warranted. Let us begin with a short review of the 
evidence for broken v^-v T symmetry. 

Following [l5j . the three-neutrino mixing matrix is 
conventionally parameterized by three planar rotations, 
analogous to the three Euler angles of classical mechan- 
ics, and a purely quantum-mechanical Dirac phase 5 [25$ : 

UpMNS = ^32(^32) Pl3(#13, S) i?l2(012) ■ (12) 

The double argument of P13 is a reminder that by con- 
vention, the quantum mechanical Dirac phase appears in 
the off-diagonal elements of the rotation in the 13-plane: 

/ cos 6x3 smO^e' 15 \ , 
Rl3 ~ {- S m9 i3 e+ iS cos0 13 J ' (13) 

The range for the angles is [0, 7r/2], while the range for 
the phase is [— it, +ir], or equivalently, cosS € [—1, +1]. 

In terms of the three angles and the single phase, 
one finds that the conditions |J7 W -| = \U T j\ for v^- 
v T symmetry require that (i) #32 = 7r/4 and (ii) 
sin(26>i 2 ) sin 0i3 cos 5 = 0. 

Choices of angles/phase which satisfy (ii) above, nec- 
essary to uphold the v^-v T symmetry, are [16J: 

• Case (a): #13 = (TBM mixing is a special sub- 
case of #13 = 0, wherein sin #12 is set to l/\/3); 

• Case (b): 9 12 = or tt/2; 

• Case (c): cos<5 = 0, i.e. 8 = ±tt/2. 

The recent spectacular evidence that #13 ~ 9° is not only 
nonzero, but many a from zero [Io| . rules out case (a). 
The value of #12 is inferred from experiment to be far 
from either zero or n/2, which rules out case (b). (In 
addition, the matter effect responsible for suppression of 
v e from the sun requires \U e i\ > \U e 2\ > 0, which also 
rules out case (b).) 



Finally, we are left to discuss case (c). There is little ex- 
perimental constraint on 5, for it occurs ( Eq. (|13p ) with 
sin #13 as a prefactor, and until this year sin #13 was con- 
sistent with zero. Recent experimental evidence [l7j sug- 
gests that #32 is not equal to the maximal-mixing value 
7r/4. Thus, case (c) is mildly disfavored by data. So it 
appears that v {i -v T symmetry is likely broken. If so, then 
P is an invertible matrix. And even if #32 were exactly 
equal 7r/4, it still remains a possibilty that v^-v T sym- 
metry is broken by a value S ^ ±7r/2. In the rest of this 
section we proceed to analyze 3-neutrino flavor propaga- 
tion in the now favored case of a (slightly) broken v^-v r 
symmetry. 

A. Three-Flavor Analysis 

The constraint W e + + W T — 1 reduces the three- 
dimensional We, W,j,, WV-space to the physical triangle 
with corners at (1,0,0), (0,1,0), and (0,0,1). Because 
the relation in Eq. © is linear in W, the extremes of w 
are obtained from the values of W at these corners. The 
result is 

max / xmn{w a } — max / min{P ae , P a/1 , P QT } ■ (14) 

(The two-component analog of this result was presented 
below Eq. ©.) 

For the theoretical expectation that v T 's are not pro- 
duced by the cosmic mechanisms, W T — and the ex- 
tremes of w are even simpler: 

max / min{w Q } = max / min{P Qe , P Q/J } , (no source v T s) 

These simply-stated results have profound meaning. For 
example, a measurement of any w a satisfying Eq. (|14p 
but not satisfying (fT5|) would establish that in fact u T 's 
are produced at cosmically-distant sources. (Note that 
there is no two-component analog of the W T — re- 
sult, because with the parameter does not distinguish 
among W T and W^.) 

Inversion of the symmetric 3x3 matrix P yields the 
symmetric inverse matrix 



Det(P) 



/ [PfjtfjtPrT Pf_LT ) {Per Pfir Pc^Ptt) (-^ejLt-^/xr PeT-Pfj.fi ) \ 
{PeePrr ~ Per) (-^er-fep, ~~ PeePfir) 



(16) 



With this P 1 matrix in hand, we may directly calculate w are measured. This matrix, along with the numerical 
the injection ratios W = P _1 w once the Earthly ratios values given below for its matrix elements, are among our 

main results. 
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As mentioned earlier, the inverse of P, namely P Q(3 1 = 
Det _1 (P) (— \) a+ P x cofactor of element P a p, does not 
exist for the TBM model, because Det(PxBM) vanishes. 
However, the sign-weighted cofactor matrix, (— 1) Q+/3 x 
cofactor of (Ptbm)c^, does exist, and its form sheds light 
on what one can expect in the realistic three-neutrino 
case for the relative values of matrix elements in P . 
From the values of the TBM matrix in Eq. (|6]), one 
can easily calculate the sign-weighted cofactor matrix, 
"Det(P)P _1 ", even for matrices with a vanishing de- 
terminant. (The quotation marks here are merely a re- 
minder that "Det(P) is just a convenient label for 
the sign- weighted cofactor matrix.) We get 



'[Det(P TB M) (Ptem)- 1 ] 



1 







sis (17) 



and we repeat for the reader that the right-hand side of 
this relation is calculated independently of Det(P). Be- 
cause the updated values of the mixing angles break the 
v^-Vt symmetry implicit in the TBM model by only a 
small amount, we may expect that with updated values 
of mixing angles, matrix element values similar to those 
of Eq. (JT7I) will result for "Det(P) P _1 " , and a texture 
similar to that in Eq. (|17[) will result for the inverse ma- 



trix P" 1 . 

Of course, with small v^-v r symmetry breaking, 
Det(P) will also be small, but as long as it is nonzero, 
the matrix P is invertible. 

It is instructive to see how the shift in Det(P) will hap- 
pen. Det(P) is unchanged when we subtract the second 
row of P from the bottom row, and then subtract the 
second column of P from the third column. After a bit 
of algebra, the result is 



Det(P) 



P, 



P, 



Wej?Aj 



p, 



Y^\U ej \ 2 A 3 



A, 



£,AA 



where Aj = \U T j\ 2 
acterizing the breaking of v 



(18) 



U^j | 2 are the three parameters char- 



two of the three Aj 



-v T symmetry [181 ]. Only 
xj are independent, since J^j Aj = 
/( ,^rj| 2 = 1 — 1 = 0. The symmetry of 
P allows the row subtraction and column subtraction in 
Det(P) to each introduce a factor of Aj, leading to a 
determinant that is of order (AjA&), very small. In the 
form of Eq. (fTg)) , the determinant is easily evaluated. The 
result is 



Det(P)= | A 



e e 



p2 \ 

efiJ 



^AjA k 



[2P e ^\U ej \ 2 \U, 



,;. ~ Pw \U ej \ 2 \U ek \ 2 - P ee \U M \ 2 |iy | 2 ) . 
(19) 



In terms of the mixing angles and Dirac phase, we note of 21] as 



that the ; 
pressions 



symmetry-breaking parameters have ex- 



A 2 
A 3 



c(20 32 ) (c 2 2 s 
c(29 32 ) (s 2 
c(20 32 )c : 



13 



2 °13 



'12 V 



13 



s(29 12 ) s{29 32 ) si3 cos 6, 
s{29 l2 ) s{29 32 ) s l3 cos 5, 
(20) 



B. Three-Flavor Numerics 

Perturbative expansions about TBM values for neu- 
trino flavor ratios and for the inference of the Dirac 
phase S have been considered previously, both before [l9| 
the measurement of 6>i3, and after 20]. None of the 
prior work considers the inversion of the propagation 
matrix. In what follows, we will not appeal to a per- 
turbative expansion; rather, we will use central values 
and errors from the most recent direct fits of the PMNS 
matrix parameters to global data. Using the notation 
[—2ct, best fit, +2cr], we summarize the global analysis 



surfl 13 C 

„:„2 , 



%2 C 



[0.019, 0.0246, 0.030], 
[0.38,0.613/0.427,0.66], (21) 



if the neutrino mass-squared ordering displays a "nor- 
mal" hierarchy (i.e., (m 2 — to 2 .) ^ (m 2 — m\) > 0); and 
slightly different, 

sin 2 6»i3 C [0.020,0.0250,0.030], 

sin 2 e» 32 C [0.39, 0.600, 0.65], (22) 

if the hierarchy is "inverted" (i.e., (m 2 — to 2 ) 3> (to 2 . — 



to 2 ) > 0). The mass ordering toi < m 2 is fixed by the 
matter effect in the Sun needed to explain the observed 
solar ratio w e ~ 1/3. This leaves the two possible hier- 
archical orderings identified above. The claim is made 
in Ref. [21| that sin 2 9\ 3 is 10.2 a away from zero. The 
first and second "best fit" options for sin 2 9 32 in the nor- 
mal hierarchy reflect an octant ambiguity in the present 
data. We remark that sin 2 9 32 is not maximal in the best 
fit value, nor within the la error, but may be maximal 
at 2a. 
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For either hierarchy, the remaining results of the global 
fit concern sin 2 8 12 and cos 5: 

sin 2 #12 C [ 0.29, 0.320, 0.35 ] , (23) 

which shows that sin 2 #12 is many a away from either 
zero or the maximal value of 1/2. For cos S, at even 
la the entire range of [— 1,+1] is allowed. We remark 
that even though our propagation matrix elements P a /3 
have a classical explanation (given in the introduction), 
they nevertheless depend on the quantum mechanical 
"Dirac phase" parameter <5, via the CP-conserving fac- 
tor Ke(e lS ) — cos 5. We will take three typical values 
(0,+lj— 1) for the unconstrained parameter cos<5. The 
global analysis of [HI, [23[ find numbers similar to those 
used here (21|. We note that for the choice cos<5 = 0, the 
deviation of #32 from maximality (tt/4) is the sole source 



of v ll -v T symmetry breaking. (Accordingly, Det(P) will 
be very small for the choice cos 5 — 0.) 

Experiments necessarily determine the above parame- 
ters in combinations. Accordingly, the parameter errors 
quoted above are correlated. We will take these errors 
as uncorrelated, since the alternative requires an inde- 
pendent global fit for each change of any parameter [24j . 
Treating the errors as uncorrelated is conservative in that 
it allows a larger range of parameter values for a given 
confidence level. 

We illustrate the result of the global best fit for the 
normal hierarchy case, with 632 < t/4. The three entries 
per matrix correspond to three assumed values for the 
Dirac CP phase, with ordering cos 6 = 0, +1,— 1. We 
find: 



0.663 0.312 0.0246 \ / 0.663 0.312 0.0246 \ / 0.663 0.312 0.0246 

\U aj \ A = ( 0.191 0.393 0.416 , 0.263 0.321 0.416 , 0.118 0.465 0.416 ] ; (24) 

0.146 0.295 0.559 / \ 0.0738 0.368 0.559 / \ 0.219 0.222 0.559 

and the symmetric matrices 

0.538 0.259 0.203 \ / 0.538 0.285 0.177 \ / 0.538 0.234 0.228 \ 
P = I • • • 0.364 0.377 , ■ ■ • 0.345 0.370 , • • • 0.404 0.362 ; (25) 
0.421 / \ 0.453 J \ 0.410 / 

0.0115 -0.0327 0.0238 \ / 0.0195 -0.0633 0.0441 \ / 0.0344 -0.0132 -0.00750 \ 
Det(P)P" 1 = ( ••• 0.185 -0.150 , ••• 0.212 -0.149 , ••• 0.168 -0.141 (26) 
0.129 / \ 0.105 / \ 0.163 / 

4.59 -13.1 9.52 \ / 66.9 -217 151 \ / 2.51 -0.961 -0.548 \ 

••• 74.2 -60.1 , ■•■ 728 -510 , ■•• 12.3 -10.3 ; (27) 
51.5 / \ 359 / \ 11.9 / 



where Det(P) = (2.50, 0.291, 13.7) x 10~ 3 . (28) 

The \U\ 2 and P matrices may be compared to the corre- 
sponding TBM matrices, given earlier in Eqs. ([5]) and ([6]). 
The determinant of Ptbm vanishes, and so there is no 
P^bm t° which one may compare. On the other hand, 
the matrix Det(P) P _1 has the form of the analogous 
TBM matrix, given in Eq. ([T7|) . 

A visual comparison between the TBM and Nature's 



choices is given in Fig. 1. The constraint w e + +w T — 
1 reduces w-space to a triangle with corners at (1,0,0), 
(0,1,0), and (0,0,1). We orient the triangle with w e at the 
apex. Then, v^-v-,- symmetry with its = w T defines 
a vertical line through the center of the triangle. The 
horizontal distance of the point w from the center line 
provides a kind of measure of v^-v T symmetry breaking. 



The left triangle plots (lower-case letters) show the en- 
tire w-parameter space, for three different flavor-injection 
models for normal (NH) and inverted (IH) neutrino mass 
hierarchies. In descending order, the plots are: 

• (a) NH with pion chain injection, 



• (b) NH with quenched muon (incomplete pion 
chain) injection, 

• (c) NH with /3-beam injection, and 

• (d) IH with /3-beam injection. 
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FIG. 1: Triangle plots of (left) entire -uJ-parameter 
space, (right) un-normalized blow-up of left panel 
parameter region. See the main text for an explanation 
of symbols. 



Each plot shows the flavor values for TBM (large 
solid dot) and nine updated sets of fitted mixing angles: 
square, star, and open circle correspond to Dirac phase 
S = 7r/2,0,7r, respectively(i.e. cos<5 = (0, 1, -1)). Larger 
symbols correspond to best fit values for mixing angles 
and smaller symbols to ±2er values. Combinations of hi- 



erarchies and injection models not shown would appear 
very similar to one of the four plots on display. Right pan- 
els (capital letters A,B,C,D) show un-normalized blow-up 
of left panel parameter region containing the nine predic- 
tions and TBM value. v^-v r symmetry predicts a value 
on the vertical line through the TBM dot; deviations 
from v^-v r symmetry are evident. 

C. Examples of Three-Flavor Phenomenology 

Another use of the numerical results is to input the 
P aj 3 values from Eq. fl2S} into Eqs. (JT3J) and fI5|). Here 
we learn, for example, that w e is bounded, in the NH 
model with #32 < t/4, by a maximum of 0.538, and by a 
minimum of (0.177, 0.203, 0.228), for cos<5 = (0,+l,-l). 
When W T is set to zero to conform with theoretical preju- 
dice, then the maximum is not affected but the minimum 
rises to (0.285, 0.259, 0.234), respectively. 

Although the main use of the inverse propagation ma- 
trix which we have constructed is implementation of 
Eq. Q to infer neutrino flavor ratios at the cosmic 
sources, here we present yet another use for the ma- 
trix P _1 . Compelling theoretical arguments from par- 
ticle physics tell us that v T production at the source is 
very suppressed, due to the heavy r mass. Thus we ex- 
pect W T ~ 0. This expectation can be easily checked. 
From Eq. (Q} we have 

= W T = P~ l w e + P^w^ + P- T V , (29) 

with the elements of P _1 given analytically in Eq. (|16p . 
Multiplied by Det(P), this result is 

= (P e pPpr-P*rPw)*»e ( 30 ) 
(Per Pen PeePfir) ^/i (PeeP^n ^*eu) ^ T ' 

subjecting to the normalization w e + Wn + w T = 1. Nu- 
merical values for the parenthetical expressions for the 
best fit of the normal hierarchy with #32 in the first octant 
are given in the final columns of Eq. (|26p . Any observed 
violation of this result would implicate v T production at 
the sources. 

As a final demonstration of the utility of the new, full 
flavor evolution matrix, we consider the dependence of 
the flavor ratio W^/We at the sources on the flavor ra- 
tios observed at Earth. Here, we embrace theoretical 
prejudice and assume that i/ T 's are not produced at the 
sources, i.e. W T — 0. We use Eq. © to derive the three 
w a {W e ,W^), then sum w e and w T to get w s h, from the 
Earthly ratio w M (W e , W^/wshiWe, W M ), and invert to 
get 

W» = Pen - {Pee + Per) 

W e (P + P ) ( ^) - P 

This equation generalizes Eq. (TTTI) to the condition of bro- 
ken v^-Vt symmetry. It is independent of any injection 
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model. Therefore, it allows us to infer the flavor ratio of 
cosmically distant sources from the track-to-shower ratio 
observed at Earth, and thereby discriminate among in- 
jection models. Values for the injection ratio expected 
from the most popular source models are ~ (2,oo,0) 
for the pion decay-chain, quenched muon, and /3-beam 
models, respectively. These values for the injection ratio 
are quite different, and so discrimination among popular 
models via Eq. (|3Tj) should be straightforward. 

IV. DISCUSSIONS AND CONCLUSIONS 

Experimental inference of flavor ratios involves some 
uncertainty. We have neglected these uncertainties in this 
more theoretical paper. For example, neutrino neutral- 
current (NC) interactions, the same for all three flavors, 
contribute to shower events. The ratio of the neutral 
to charged current cross-sections is known, so the NC 
contribution can be accounted for in a data sample. As 
another example, the experimental efficiencies for mea- 
suring shower events and muon-track events are different. 
Again, these can be accounted for in a data sample. 

What we have shown is that neutrino flavor physics 
is rich in its information content, and therefore worth 
pursuing. Neutrino flavor physics offers another window 
into the dynamics of the most distant, most energetic 



objects in the Universe. 

We end with a summary of the main results discussed 
in this paper. We have: 

• shown the movement of Earthly measured flavor 
ratios away from the v^-Vt symmetric value of 
the previously viable TriBiMaximal model, for the 
three most popular cosmic-source flavor models; 

• derived the inverse flavor propagation matrix 
which allows one to infer flavor ratios injected at 
cosmically-distant sources from the ratios observed 
here on Earth; 

• presented a relation among flavor ratios at Earth 
that determines whether i/ T 's are injected at the 
source; 

• derived a general formula for the v^jv^ injection 
flavor ratio at the source in terms of the observable 
ratio of track-to-shower events at Earth. 
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